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REAL-TIME SIMULATION OF FlOO-PW-lOO TURBOFAN ENGINE 
USING THE HYBRID COMPUTER 
by John R. Szuch and Kurt Seldner 
Lewis Research Center 

SUMMARY 

This report describes the development of a real-time hybrid computer simulation of 
the FlOO-PW-lOO augmented turbofan. The simulation is intended to support controls 
research programs involving that engine. The FlOO-PW-lOO real-time simulation was 
modeled after a similar simulation of the TF30-P-3 turbofan engine. The FlOO-PW-lOO 
real-time simulation has both wide-range steady-state and transient computing capabil- 
ities. Factors such as fluid momentum, mass and energy storage, and rotor inertias 
are included. The digital portion of the hybrid computer was used primarily to perform 
the bivariate function generation associated with modeling the performance of the en- 
gine's rotating components. The remaining calculations (including integration) were per- 
formed by the analog portion of the hybrid computer . 

Both steady-state and transient data were generated with the real-time engine sim- 
ulation. An analog simulation of the engine's hydromechanical fuel control and a digital 
simulation of the engine’s electronic supervisory control were used to control the sim- 
ulated engine over a wide range of power settings and flight conditions. Resultant 
steady -state data were compared with similar data obtained with the engine manufac- 
turer’s baseline digital simulation of the engine and controls. The results of the com- 
parison show that the real-time simulation results do adequately match the baseline 
digital steady -state results over the flight envelope. The responses of both the real-time 
hybrid and baseline digital simulated engines to rapid movements of the power lever were 
compared at sea -level static conditions. The results of that comparison indicate that the 
real-time simulation does adequately match the baseline digital simulation in predicting 
the transient behavior of the engine. Steady-state simulation data were compared with 
sea-level, static experimental data. The results of this comparison indicate that both 
the baseline digital and real-time hybrid simulations of the FlOO-PW-lOO engine do ad- 
equately predict the performance of the actual engine. 



INTRODUCTION 


In recent years, there has been increased interest in developing digital, electronic 
controls for airbreathing propulsion systems (ref. 1). The use of a digital computer, 
with its inherent precision and its logic and memory capabilities, will provide the pro- 
pulsion system control with more flexibility and versatility than is currently provided by 
hydromechanical controls. The increased flexibility and versatility are necessary to al- 
low engines and inlets to operate at higher levels of performance over a wider range of 
flight conditions. The use of digital controls can also reduce pilot workloads, improve 
control maintainability, and reduce control complexity and weight. 

The development of digital controls for turbojet and turbofan engines can be facil- 
itated by the use of a real-time computer simulation of the engine (ref. 2). The engine 
simulation provides a ’’test bed" for evaluating new control laws and for checking and 
"debugging" the actual control software prior to engine testing. A real-time engine 
simulation enables the control developer to evaluate the timing and sequencing within the 
digital control and to predict the effects of extended digital sampling intervals (refs. 3 
and 4) on engine performance. 

Turbojet and turbofan engines may be simulated by using analog (ref. 5), digital 
(refs. 6 to 8), or hybrid (refs. 2, 9, and 10) computers. A real-time requirement usu- 
ally limits the possible choices to the analog and hybrid computers. A purely analog 
simulation of a turbofan engine requires a large amount of computing equipment. Much 
of this equipment is needed to perform the function generation associated with modeling 
the performance of the engine's rotating components (fan, compressor (or compressors), 
and turbine (or turbines)). Large, all -analog simulations may also suffer from long set- 
up and checkout times and a lower level of accuracy. The precision, repeatability, and 
memory capabilities of a hybrid computer can be used to perform the necessary function 
generation and also to aid in the setup and checkout of the analog portion of the 
simulation. 

This report describes a real-time hybrid computer simulation of the Pratt & 

Whitney F100-PW-100 augmented turbofan engine. The simulation has both steady-state 
and transient calculation capabilities and is intended for supporting controls research 
programs involving that engine. The mathematical model of the engine was implemented 
on the Lewis Research Center's Electronic Associates Incorporated (EAI) Hybrid Com- 
puter System. This system consists of a Model 640 Digital Computer, a Model 680 
Analog Computer, and a Model 681 Analog Computer. The mathematical model was pat- 
terned after the engine manufacturer's digital (baseline) simulation of the engine. Many 
of the assumptions and techniques employed in the development of the F100-PW-100 real- 
time simulation are based on experience gained while developing a similar real-time 
simulation of the Pratt & Whitney TF30-P-3 engine (ref. 2). 
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Results are presented which compare steady -state data from the real-time simula- 
tion with baseline digital simulation results over a wide range of power settings and flight 
conditions. Limited comparisons with transient results from the baseline digital simula- 
tion and with steady -state experimental data are also made. The report also includes the 
equations which model the engine, FORTRAN listings, and analog patching diagrams. 
Significant differences between the F100-PW-100 and TF30-P-3 real-time simulations 
are pointed out. 


ENGINE DESCRIPTION 


The Pratt & Whitney F100-PW-100 engine (fig. 1) is an axial, mixed -flow, aug- 
mented, twin-spool, low -bypass -ratio turbofan. A single inlet is used for both the fan 
airflow and the engine core airflow. Airflow leaving the fan is separated into two flow 
streams: one stream passing through the engine core and the other stream passing 
through the annular fan duct. A three -stage fan is connected by a through-shaft to the 
two-stage, low-pressure turbine. A 10-stage compressor is connected by a hollow shaft 
to the two-stage, high-pressure turbine. The fan has variable, trailing -edge, inlet guide 
vanes. The inlet guide vanes are positioned by the engine's electronic supervisory con- 
trol system as a function of fan corrected speed to maintain fan stability at low speeds. 
The compressor has a variable inlet guide vane followed by two variable stator vanes. 

The compressor vanes are positioned by the engine's hydromechanical fuel control sys- 
tem as a function of compressor corrected speed. Engine airflow bleed is extracted at 
the compressor exit and discharged through the fan duct during starting. Compressor 
discharge bleed air is also used to cool the high- and low-pressure -turbine blades and to 
power the augmentor turbopump. 

The main combustor consists of an annular diffuser and a chamber with 16 fuel noz- 


zles. The hydromechanical fuel control system meters fuel to the main combustor as a 
function of the power lever angle PLA, the compressor speed N H , the fan discharge 
total temperature T^g, and the compressor discharge static pressure P g g. (All sym- 
bols are defined in appendix A. Numerical subscripts refer to locations in the engine 
(fig. 1). ) 

The engine core and fan duct streams combine in an augmentor and are discharged 
through a variable convergent -divergent nozzle. The augmentor consists of a diffuser 
section and five concentric fuel manifolds (zones). The hydromechanical augmentor fuel 


control meters fuel to the five zones as a function of PLA, T^, and P g g. The nozzle 
area is controlled so as to maintain a desired engine airflow during augmented operation. 
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ENGINE MODEL 


The mathematical model describing the performance of the F100-PW-100 engine is 
patterned after the engine manufacturer's digital simulation (CCD 1015) of the engine. 
Wide-range, overall performance maps of the engine's rotating components are utilized 
so as to provide wide-range, steady-state accuracy, Factors such as fluid momentum, 
mass and energy storage, and rotor inertias are included so as to provide transient ca- 
pability. However, in order to satisfy the real-time requirement and to minimize the 
amount of required analog equipment, nearly all fluid properties such as specific heats 
and gas constants are assumed to be constant. 

Figure 2 contains a computational flow diagram of the F100-PW-100 real-time sim- 
ulation. It can be seen that the mathematical model consists of a number of individual 
elements and their related volumes (inlet, fan, compressor, etc. ) each of which re- 
quires a number of input variables and generates one or more output variables. How- 
ever, unlike a digital simulation, all calculations in the real-time hybrid simulation can 
be considered to be performed in parallel. Descriptions of the individual element 
models are provided in appendix B. Appendix B also contains a complete list of equa- 
tions defining the F100-PW-100 real-time simulation model. Differences between the 
F100-PW-100 and TF30-P-3 (ref. 2) models are pointed out. 


HYBRID SIMULATION 

The equations describing the F100-PW-100 mathematical model were implemented 
on the Lewis Research Center’s Hybrid Computing System. This system consists of an 
EAI Model 640 Digital Computer, a Model 680 Analog Computer, and a Model 681 Analog 
Computer. The split of the computation load between the digital and analog computers 
was similar to that employed in the TF30-P-3 simulation (ref. 2). That is, the digital 
computer was used primarily to perform the bivariate function generation associated 
with modeling the performance of the engine’s rotating components. In addition to this 
function generation, the digital computer also performed the calculations required to 
compute the engine thrust. In order to minimize the core storage requirement and the 
digital update time while allowing the use of FORTRAN programming, scaled -fraction 
(ref. 2) variables and arithmetic routines were used in the main digital computational 
loop. Scaled -fraction variables may not exceed unity and are represented by a single 16- 
bit word on the EAI 640 Digital Computer. The need to limit the digital update time to a 
small value (~5 msec) to guarantee stability when operating the simulation in real time 
prevented a fuller use of the digital computing capability. Appendix C contains a 
FORTRAN listing of the digital portion of the F100-PW-100 real-time simulation. Ref- 
erence 2 includes a discussion of the digital function generation programs MAP 2 and 
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MAP2L, which are used in both the TF30-P-3 and F100-PW-100 simulations. The dig- 
ital portion of the F100-PW-100 real-time simulation required 11 741 words of core 
storage and had an update time of 6. 37 milliseconds with constant values for the input 
variables. 

The remaining calculations were performed on the analog computers. The analog 
computers performed all the operations characteristics of analog computers (i.e. , sum- 
ming, integration with respect to time, multiplication, attenuation, univariate function 
generation, etc. ). The use of peripheral equipment such as x-y plotters and strip-chart 
recorders allowed continuous monitoring by the user of computed variables. The split of 
the computation load between the two analog consoles was based on the computing equip- 
ment limitations of each console. For example, the full complement of 24 multipliers on 
the 680 Analog Computer and 30 multipliers on the 681 Analog Computer was utilized. In 
addition, the full complement of six digital -analog multipliers (DAM’s) on the 680 Analog 
Computer was used. Appendix D contains the analog patching diagrams and potentiom- 
eter settings for the F100-PW-100 real-time simulation. 


RESULTS AND DISCUSSION 

The usefulness of a simulation such as the F100-PW-100 real-time simulation de- 
pends on its ability to accurately represent the physical system over the desired range of 
operation. The F100-PW-100 real-time simulation should match the steady -state and 
transient engine performance for power settings from idle to full -augmentation at alti- 
tudes and flight speeds within the envelope shown in figure 3. Figure 3 also shows the 
flight conditions selected for comparison of hybrid and baseline digital data. 

The following sections compare both steady-state and transient data obtained with 
the F100-PW-100 real-time simulation with the corresponding baseline digital data. The 
validity of both the hybrid and baseline digital simulations is demonstrated by comparing 
steady-state simulation data with available experimental data. 


Steady -State Simulation Results 

Comparison with baseline digital data . - The F100-PW-100 real-time simulation 
was operated in an open -loop fashion at each of the flight conditions shown in figure 3. 
Open-loop operation was selected to allow any engine -related errors to be isolated from 
control -related errors. For all open-loop runs, the five input variables to the simula- 
tion (Wp 4> w p 7, A„, GVIPOS, and HVSPOS) were set to the values obtained from the 
baseline ’digital simulation at the flight condition and power setting of interest. This pro- 
vided a basis of comparison between hybrid and baseline digital data. The engine 
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variables selected for comparison were compressor speed, fan speed, main -combustor 
pressure, high-pressure-turbine inlet temperature, fan airflow, fan -tip pressure ratio, 
and nozzle inlet temperature. Agreement between hybrid and baseline digital values for 
all these variables implies agreement for almost all other engine variables. For con- 
venience, engine variable data were plotted against power lever angle. 

Figures 4 to 10 show the effects of flight condition on the match of hybrid and base- 
line digital steady -state data. Since the ’’design" point for the hybrid simulation was the 
sea-level, static, military -power (PLA = 73°) condition, the agreement between hybrid 
and baseline digital data at this condition was expected. Therefore, the validity of the 
mathematical model used in the hybrid simulation depends on the match of hybrid and 
baseline digital data at part-power settings and/or non-sea-level, static conditions. 

Figures 4 and 5 compare the hybrid and baseline digital simulation values of com- 
pressor and fan speeds, respectively. Good agreement between rotor speeds at part- 
power settings for all flight conditions was observed. The observed disagreement be- 
tween fan speeds during non-sea-level, static augmentation (PLA > 80°) may be attrib- 
uted to differences in the hybrid and baseline digital models of the nozzle. The assump- 
tion of constant, nozzle inlet, specific-heat ratio in the hybrid simulation would contrib- 
ute to the errors in fan speed, which amount to about 4 percent of the design value. 

Figures 6 and 7 compare the hybrid and baseline digital simulation values of main- 
combustor pressure and high -pressure -turbine inlet temperature, respectively. The ob- 
served agreement in pressure and temperature over the range of power settings and 
flight conditions substantiates the steady -state models of the fan hub, compressor, and 
main combustor used in the hybrid simulation. 

Figures 8 and 9 compare the hybrid and baseline digital simulation values of fan air- 
flow and fan-tip pressure ratio, respectively. Some discrepancies in both airflow and 
pressure ratio were observed at the altitude conditions. These discrepancies may be due 
to a combination of (1) no Reynolds number effects on fan performance in the hybrid sim- 
ulation and (2) the aforementioned differences in the nozzle models at non -sea -level, 
static conditions. The observed discrepancies in fan airflow and fan -tip pressure ratio 
at high altitudes amount to approximately 3 percent of the corresponding design values. 

Figure 10 compares the hybrid and baseline digital simulation values of the nozzle 
inlet temperature. Reasonable agreement was observed over the range of power settings 
and flight conditions. The observed differences in nozzle inlet temperature at the highest 
altitude conditions may be attributed to the assumption of constant specific heats in the 
hybrid computer model of the augmentor. 

Comparison with experimental data. - In general, figures 4 to 10 indicate that the 
hybrid simulation does adequately match the baseline digital simulation in representing 
the steady-state behavior of the F100-PW-100 engine. The ability of both the hybrid and 
digital simulations to represent the actual engine performance, however, can only be 
demonstrated by comparing simulation data with experimental data. Steady -state , 
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experimental data obtained at sea-level, static conditions were available. Since the non- 
augmented and augmented tests were run at different ambient temperatures, data from 
those tests were considered separately. The hybrid simulation was run with the appro- 
priate inlet conditions to match the experimental tests. The controlled inputs to the sim- 
ulation were provided by an analog simulation of the hydromechanical fuel control and a 
digital simulation of the electronic supervisory control. Because of a discrepancy be- 
tween the experimental and simulated values of the power lever angle, comparisons of 
hybrid, baseline digital, and experimental data were made on the basis of equivalent 
main -combustor fuel flows (for nonaugmented tests) and augmentor fuel flows (for aug- 
mented tests) . 

Figure 11 compares simulation and experimental data for the nonaugmented tests. 

The excellent agreement in all engine variables indicates (1) that the baseline digital sim- 
ulation does represent the actual engine sea -level static performance over the entire 
range of nonaugmented operation and (2) that the hybrid simulation does adequately match 
the baseline digital simulation. 

Figure 12 compares simulation and experimental data for the augmented tests. 

Since most engine variables are held at their military-power values by the control during 
augmentation, the only significant variables to consider during augmentation are the noz- 
zle inlet temperature and the nozzle area. Figure 12 indicates that the simulations do 
predict the correct nozzle inlet temperature resulting from the augmentor fuel flow. 
However, figure 12(b) shows that the actual engine control resulted in a larger nozzle 
area during augmentation than predicted by either simulation. Figure 12 indicates that 
the actual nozzle may have a lower flow coefficient than predicted by the baseline digital 
simulation. 


Transient Simulation Results 

The preceding sections have demonstrated the ability of the F100-PW-100 real-time 
hybrid simulation to predict the steady -state performance of the actual engine. If the 
simulation is to serve as a useful tool in developing digital control systems for the en- 
gine, it must also be capable of predicting transient engine performance. A comparison 
of hybrid and baseline digital transient data was limited to sea-level, static conditions 
because additional baseline digital data were not available. 

The hybrid simulation was operated with the five control inputs scheduled as func- 
tions of time to match the baseline digital values for a power lever ramp from idle to 
military power in 0. 5 second. As in the case of the steady -state comparison of simula- 
tion results, open-loop operation was selected for the transient comparison to isolate any 
simulation errors from errors associated with the control implementation. 
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Figure 13(a) shows a comparison of the hybrid and baseline digital simulation re- 
sponses of the main -combustor pressure. The hybrid simulation response was more 
rapid than the digital response. Similar results were obtained with the real-time hybrid 
simulation of the TF30-P-3 engine (ref. 2). In that instance, the gain of the main- 
combustor integrator was decreased by a factor of 25 to match baseline digital data. 

This decrease, in effect, accounted for heat -transfer effects in the main combustor, 
which were not included in the hybrid model. Figures 13(b) and (c) show the effect of re- 
ducing this gain on the hybrid simulation response of the F100-PW-100 main -combustor 
pressure. A reduction factor of 20 was selected for subsequent transient tests. 

Figure 14 shows a comparison of the responses of hybrid and baseline digital sim- 
ulation engine variables to the simulated idle -to -military power lever "slam. ” As pre- 
viously noted, the main -combustor pressure gain was reduced by a factor of 20 for this 
test. Also, the scheduled values of exhaust nozzle area correspond here to an "idle- 
area -reset on" condition. That is, the exhaust nozzle is near its maximum opening at 
idle and is scheduled closed by the hydromechanical control as a function of a rate- 
limited power lever angle as the engine accelerates. Figure 14 shows excellent agree- 
ment in all engine variables. 

Figure 15 shows a comparison of hybrid and baseline digital simulation variables for 
a simulated power lever cycle from military to maximum augmentation to military. As 
in the case of the idle -to -military power lever transient, the scheduled values for the 
control inputs were determined from the baseline digital simulation. The function of the 
exhaust nozzle area control during augmentation is to maintain a constant flow impedance 
to the main engine while the main fuel control modulates fuel flow to the main combustor 
to maintain constant compressor speed. The results shown in figure 15 indicate that the 
scheduled values for the control inputs do result in a match of baseline digital and hybrid 
simulation values for rotor speeds, main -combustor pressure and temperature, and noz- 
zle inlet temperature. However, the suspected shortcomings in the exhaust nozzle model 
did result in steady -state errors in the hybrid -simulated exhaust nozzle inlet pressure 
during maximum augmentation. No appreciable errors in the transient characteristics 
of these variables can be observed in this figure. 


SUMMARY OF RESULTS 

This report describes a real-time, hybrid computer simulation of the F100-PW-100 
augmented turbofan. While the F100-PW-100 real-time simulation was modeled after the 
TF30-P-3 simulation, a number of modifications have been incorporated in an attempt to 
extend the range of accuracy. For example, temperature -sensitive specific heats are 
used in the calculations of fan and compressor torques. Steady-state simulation data are 
presented for a full range of power settings over the entire flight envelope. The 
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F100-PW-100 real-time simulation is shown to adequately match the baseline digital 
steady -state results over the entire flight envelope. Transient data are presented for 
sea-level, static conditions and are also shown to match digital results. Both hybrid 
and baseline digital steady -state data are compared with sea-level, static experimental 
data to show that the basic mathematical models used in the simulations are valid. 
These results indicate that the real-time simulation requirement can be satisfied while 
maintaining a sufficient degree of accuracy over a wide range of power settings and 
flight conditions. 

With the simulation structure described in this report the digital computer update 
time was 6.37 milliseconds and resulted in stable, real-time operation. The use of a 
faster digital computer than the EAI Model 640 would, of course, allow more computa- 
tions to be done in the digital computer during the 6- to 7 -millisecond duration. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, March 11, 1975, 

505-04. 
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APPENDIX A 


SYMBOLS 
2 2 

A cross-sectional area, cm (in. ) 

nozzle flow coefficient 

c specific heat at constant pressure, J/kg-K (Btu/lbm-°R) 

r 

F thrust, N (lbf) 

FN7 nozzle flow function 

fj functional relation, i = 1 to 11 

f/a local fuel -air ratio 

GVIPOS fan inlet guide vane position, deg 

g c gravitational conversion factor, 100 cm-kg/N-sec 2 (386. 3 Ibm-in. /lbf-sec 2 ) 

HVF heating value of fuel, J/kg (Btu/lbm) 

HVSPOS compressor stator vane position, deg 
Ah turbine enthalpy drop, JAg (Btu/lbm) 

hp turbine map enthalpy drop parameter, JAg-K^-rpm (Btu/lbm ^R 1 / 2 -rpm) 

I polar moment of inertia, N-cm -sec (ft -lbf-sec ) 

J mechanical equivalent of heat, 100 N-cm/J (778.3 ft-lbf/Btu) 

K AB augmentor pressure loss coefficient, 

N 2 -sec 2 Ag^ -cm^ -K (lbf 2 -sec 2 /lbm 2 -in. ^ -°R) 

Kg main -combustor pressure loss coefficient, 

N 2 -sec 2 Ag^ -cm^ -K (lbf 2 -sec 2 /lbm 2 -in. ^ -°R) 

K BLWHT frac ^ ion of high -pressure -turbine cooling bleed that is performing work 

Krlwlt ^ rac ^ on low-pressure -turbine cooling bleed that is performing work 

K n nozzle flow constant, kg-K 1//2 /N-sec (lbrn^R^ 2 /lbf-sec) 

K pR5 low -pressure -turbine discharge pressure loss coefficient 

l length, cm (in. ) 

M Mach number 

N rotational speed, rpm 

2 

P total pressure, N/cm (psia) 
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P/P pressure ratio 

PLA power lever angle, deg 

o 

P g static pressure, N/cm (psia) 

Q torque, N-cm (in. -lbf) 

gas constant of air, 44. 83 N-cm/kg-K (640. 1 in. -lbf/lbm-°R) 

T total temperature, K (°R) 

T/T temperature ratio 

t time, sec 

V volume, cm (in. ) 

W stored mass, kg (lbm) 

w mass flow rate, kg/sec (lbm/sec) 

w c corrected mass flow rate, kg/sec (lbm/sec) 

w turbine map flow parameter, kg -K -cm^/N -rpm -sec (lbm-°R-in. ^/lbf-rpm-sec) 

P 

y specific -heat ratio 

6 total pressure relative to sea-level conditions 

77 efficiency 

6 total temperature relative to standard -day conditions 

t time constant, sec 

Subscripts: 

AB augmentor 

B main combustor 

BLHT high-pressure -turbine cooling bleed 
BLLT low-pressure -turbine cooling bleed 
C compressor 

D fan duct 

e nozzle exit plane 

F fuel 

FAN fan 

H high 

HT high-pressure turbine 
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I inlet 

ID fan hub (core) 

i initial conditions 

j engine station (fig. 1); j = 0, 2, 2.1, 2.2, 3, 4.1, 5, 6, 7, 8, 13, 16 

j’ entrance to volume at station j ; j = 3, 4, 4. 1, 6, 7, 13, 16 

L low 

LT low-pressure turbine 

M map 

m measured 

N nozzle 

n net 

OD fan tip (bypass) 

A 

SUB subsonic 

SUP supersonic 

TPBL turbopump bleed 
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APPENDIX B 


ENGINE MODEL 

This appendix describes the significant features of the individual elements which 
comprise the F100-PW-100 engine model shown in figure 2. Emphasis is placed on not- 
ing the differences between the F100-PW-100 model and the TF30-P-3 model described 
in reference 2. 


Inlet 


The mathematical treatment of the inlet in the F100-PW-100 real-time simulation is 
identical to that used in reference 2. A steady -state representation of a typical inlet re- 
covery is used to provide the proper fan inlet conditions (Pg.Tg) for a specific flight con- 
dition (Pq,Tq,Mq). The following equations define the inlet model in the F100-PW-100 
real-time simulation: 


(?) 


|^-1.0 + 0.2M» 


(Bl) 


(DC 


(B2) 


77j = 1. 0 if M Q < 1. 0 


= 1.0 - 0.075(M Q - 1.0) 


1.35 


(D 


T 2 = li) T 0 


otherwise 


(B3) 

(B4) 


P 2 %) Vi 


(B5) 
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Fan 


The mathematical treatment of the fan in the F100-PW-100 real-time simulation is 
quite similar to that discussed in reference 2. That is, overall performance maps 
(rather than individual stage data) are used, with separate performance maps for the tip 
(bypass) and hub (core) sections of the fan. Figure 16 shows the fan performance maps 
with the inlet guide vanes at their axial position. Since the fan’s performance is affected 
by the inlet guide vanes, the data shown in figure 16 must be adjusted for nonaxial inlet 
guide vane positions. Based on data from the baseline digital simulation, the effect of 
the variable inlet guide vanes was accounted for by multiplying the map value of the cor- 
rected airflow (fig. 16(a)) by a linear function of the inlet guide vane position to obtain 
the actual fan corrected airflow. The effect of the inlet guide vane position on the fan- 
hub performance (fig. 16(b)) was assumed to be negligible. 

The baseline digital data also indicated that the fan -tip and fan -hub temperature 
ratios could be computed from piecewise linear functions of the corresponding pressure 
ratios. Three segments were required to fit each of the temperature -ratio - pressure- 
ratio curves. The use of the linear functions eliminates the need for the bivariate func- 
tion generation and exponentiation associated with using isentropic efficiency maps. The 
fan efficiency (hence, temperature ratios) was assumed to be unaffected by changes in 
the inlet guide vane position. The following equations define the fan model in the F100- 
PW-100 real-time simulation: 


p n n t 

(w ) = f A , -= 

c FAN.M 1 \P 2 Jf 2 


= f, 


rP 13 N L 


'FAN, ID 


2 1 P n ’ Ale, 


W 2 = (w c ) 


&j(l. 0143 + 0. 002653 GVIPOS) 


FAN.M 


P 2. 1 " P 2.2 


'FAN, ID 


(B6) 


(B7) 


(B8) 


(B9) 
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Po 


T ' = 0. 1833 Hi + 0. 9029 if — li>1.85 


P, 


P„ P13 

= 0.2533 — + 0.7733 if 1.4 < — <1.85 


0. 3200 Hi + 0. 6800 otherwise 
P 2 


(BIO) 


T 13* “ 


(?) 


FAN.OD 
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+ 0.7021 otherwise 
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(B13) 


FAN. ID 


The primed station designations (e. g. , 13 ) refer to the entrance to the volume. 


Compressor 

The F100-PW-100 real-time compressor model differs from the compressor model 
described in reference 2 in a number of respects. While both models use overall per- 
formance maps, the F100-PW-100 compressor model uses pressure ratio (rather than 
corrected airflow) as an input variable (fig. 2). This eliminates the need for an external 
calculation of the airflow. Figure 17 shows the F100-PW-100 compressor performance 
with the stator vanes at their axial position. As in the case of the fan, the performance 
of the compressor was affected by the variable stator vanes. Based on baseline digital 
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data, the effect of the variable stator vanes was accounted for by adjusting the map 
value of the corrected airflow with a bivariate function of the stator vane position and 
compressor corrected speed. Figure 18 shows the effect of variable stator vane position 
on the compressor performance map. 

The F100-PW-100 compressor temperature ratio was fit by a two-segment linear 
function of the compressor pressure ratio. The compressor efficiency (hence, temper- 
ature ratio) was assumed to be unaffected by changes in the stator vane position. The 
following equations define the F100-PW-100 real-time compressor model: 


( w c } = f 3 
c , m 2 


P 3 N H 


.2 


W 2.2 = ^c ] 


J 2.2 


C , m 


N H 

1. + f 4 ( n , HVSPOS 




s )j 


y»2” 


(B14) 


(B15) 



0.0702 



1.3814 



5.75 


= 0. 1018 


2.2 


T 3’ = 


+ 1. 1996 otherwise 


a 


2.2 


(B 16 ) 

(B17) 


Engine Bleeds 

The extraction of turbine cooling and augmentor turbopump bleeds from the com- 
pressor discharge is accounted for in the F100-PW-100 real-time simulation. All 
bleed flows were assumed to be proportional to the compressor inlet airflow. For con- 
venience, the flow loss due to compressor seal leakage was lumped with the low- 
pressure -turbine cooling bleed. The following equations define the engine bleed model 
in the F100-PW-100 real-time simulation: 

W BLHT = ^2 2 (B18) 
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^BLLT' 0 - 03828 ^^ (B19) 

w T pbl = 0.00353 Wg ^ (B20) 


Turbines 

The mathematical models describing the performance of the F100-PW-100 turbines 
are identical to those described in reference 2. The overall performance of each turbine 
is described by a pair of bivariate functions. Figures 19 and 20 are performance maps 
showing the turbine flow and enthalpy drop parameters for the high- and low-pressure 
turbines, respectively. The cooling bleed for each turbine was assumed to reenter the 
flow stream at the discharge of the turbine, although a proportion of each bleed was as- 
sumed to perform turbine work. The following equations define the turbine models in the 
F100-PW-100 real-time simulation: 


(wJ =f E 


fp 4.1 N H 


P HT 5 \ P, ’ */T 


w 4 = (W p ) 


P 4 N H 
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(hp) H T - f 6 


' P 4.1 N H 


P 4 V t 4; 


(^h)RT - ^P^HT V P 4 


(wT) 


4. 1' = w 4 


c p,4 T 4 


+ W BLHT 

c p,3 T 3 n-BLWHT^HT 

_ C P,4. 1 

_ 

c p,4.1_ 

_ C P ,4. 1 C P,4. 1 


p s N t 

(w ) =f 7 [ , - 

Plt Vp 4 .i 


*4. 1 = ( V LT P 4- 1 


N t 


4.1 


(B21) 

(B22) 

(B23) 

(B24) 

(B25) 

(B26) 

(B27) 
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(B29) 
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c p > 3 T 3 ^LWLT^^LT 


L c p ,6 


P> 6 


(B30) 


Table I lists the design parameters (specific heats, k blwhT’ K BLWLT’ etc - ^ used in 
the real-time simulation of the F100-PW-100 engine. 


Main Combustor and Augmentor 


The only significant difference between the F100-PW-100 real-time main-combustor 
model and that described in reference 2 is the elimination of the high-frequency flow dy- 
namics that were included in the TF30-P-3 simulation. The elimination of the main- 
combustor flow dynamics resulted in the need for a direct calculation of the main- 
combustor airflow from the instantaneous pressure drop across the main combustor. A 
simple energy balance was used to compute the steady-state temperature rise in the 
F100-PW-100 main combustor. However, in order to match baseline digital data over 
the full range of power settings, it was necessary to include a temperature-sensitive 
energy release from the injected fuel. The combined effects of temperature on the main- 
combustor efficiency and the heating value of the fuel were achieved with a linear function 
of the main -combustor temperature. The following equations define the main -combustor 
model in the F100-PW-100 real-time simulation: 


*3 * X- 


'P 3 (P 3 - P 4 ) 


K B T 3 


V K B V 3 
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(wT) 4 » = -EilwgTg + JJg w 
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(B32) 


-P,4 


“P,4 


HVF 

tj b £Ll£ = 49 602 - 7. 3376 T 4 
C P,4 


(B33) 
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The most significant difference between the F100-PW-100 and TF30-P-3 real-time 
augmentor models is the mixing of the core and duct flows at the augmentor entrance in 
the F100-PW-100 real-time simulation. While both real-time simulations assume a 
static pressure balance at this point in the engine, the TF30-P-3 real-time hybrid model 
(ref. 2) assumes no flow mixing and treats the core and duct sections of the augmentor 
and the nozzle as separate elements. 

In order to match F100-PW-100 baseline digital data, the effective augmentor effi- 
ciency was fit by a piecewise linear function of the fuel -air ratio. Figure 21 shows the 
augmentor efficiency required to match the baseline digital data. The following equations 
define the augmentor model in the F100-PW-100 real-time simulation: 


P 16 “ P 6 

(B34) 

P 5 = K PR5 P 6 

(B35) 

r _ p y ,: ; 2 Te - p K AB V 6 w 6 
’ 6 - K AB w 6 - P 6 RaW(s 

(B36) 

c n fi • HVF • 

(wT) 7 t~ ’ W 6 T 6 +7J AB w p,7 

(B37) 

C P,7 C p,7 


,AB=f9 (a), 

(B38) 


/f\ ^F,7 

W 7 *6"*F,4 


(B39) 


Duct 

The F100-PW-100 real-time duct model is essentially the same as the model de- 
scribed in reference 2. However, the steady -state duct pressure drop has been function- 
ally related to the duct airflow (corrected to the duct discharge conditions) in the F100- 
PW-100 simulation. Figure 22 shows the duct pressure loss characteristic. As in the 
case of the TF30-P-3 simulation, the temperature rise (or drop) in the F100-PW-100 
duct has been neglected. The following equations define the F100-PW-100 real-time duct 
model: 
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c p, 16 T 16 = c p, 13 T 13 


(B41) 


Exhaust Nozzle 

As in the case of the augmentor models, the F100-PW-100 and TF30-P-3 exhaust 
nozzle models differ with regard to the mixing of core and duct flow streams. Because 
of the assumed mixing of streams in the F100-PW-100 augmentor, only one nozzle model 
is required. 

Because of the complexity involved in computing both subsonic and supersonic flows 
in a convergent -divergent nozzle, certain approximations and correction factors were 
employed in calculating the F100-PW-100 nozzle flow. The flow was calculated as if the 
nozzle were a convergent -only design. The resulting error in flow was compensated for 
by the use of a variable nozzle flow coefficient. In order to match baseline digital data 
(nonaugmented) , a two -segment linear function of the nozzle pressure ratio was em- 
ployed. In addition, a temperature -sensitive reduction in the effective nozzle throat 
area was used to match baseline data for augmented operation. Figure 23 shows the re- 
quired nozzle area reduction. This area reduction was fit by four linear functions of the 
nozzle inlet temperature. 

Since it was desirable to have a nozzle thrust calculation in the F100-PW-100 real- 
time simulation, it was necessary to account for the effects of the divergent nozzle on 
the exhaust -plane static pressure. Based on available digital, data, the expansion ratio 
of the nozzle was fit by a quadratic function of the variable nozzle throat area. Isen- 
tropic flow tables (ref. 11) were implemented by quadratic functions of the expansion 
ratio to allow a determination of the static pressure at the nozzle exit for both subsonic 
and supersonic flow in the nozzle. Ram drag effects on the engine net thrust were also 
included in the F100-PW-100 real-time simulation. The following equations define the 
F100-PW-100 exhaust nozzle model: 

(B42) 
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w 7 =K N P 7 (FN7)C d;7 

f n =0.0 if T 7 < 787.8 

= 0.28719 T ? - 226.2 if 787. 8 < T ? < 888. 9 

= 0.09784 T ? - 57.93 if 888. 9 < T ? < 1165 

= 0. 20032 T ? - 177. 5 if 1165 < T ? < 1496 

= 0. 30867 T ? - 339. 5 if T 7 > 1496 

Ao o 

— = 0.52152 + 3.16282 A ? -2.9508 A~ 
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Engine Dynamics 

The previous sections describe the steady -state performance of the various com- 
ponents which comprise the F100-PW-100 engine. However, a transient engine simula- 
tion must also account for the effects of fluid compressibility, fluid momentum, energy 
storage, and rotor inertias on the transient performance of the engine. 

Intercomponent volumes were assumed at engine locations where either (1) gas dy- 
namics were considered to be important or (2) gas dynamics were required to avoid the 
need for iterative solution of equations. In these volumes, the storage of mass and en- 
ergy occurs. The dynamic forms of the continuity, energy, and state equations were 
solved for the stored mass, temperature, and pressure in each volume. When mixing of 
gases was not involved, a simple first-order lag form of the energy equation was used. 

The effects of fluid momentum on the transient behavior of the F100-PW-100 engine 
were considered only in the fan duct and augmentor tailpipe. The contribution of flow dy- 
namics in the compressor, main combustor, and turbines was assumed to be primarily 
high frequency (>10 Hz) in nature and was consequently ignored. 

The most significant factors in determining the transient behavior of a turbofan en- 
gine are the rotor moments of inertia. Rotor speeds were computed from the dynamic 
form of the angular momentum equation. 

A comparison of hybrid and baseline digital data at low power settings indicated that 
the assumption of constant fluid properties in the hybrid simulation can lead to simulation 
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errors. Similar results were noted in reference 2. Therefore, the calculation of the 
F100-PW-100 compressor torque requirement was expanded to include a temperature- 
sensitive specific heat at the compressor discharge. A linear function of temperature 
was found to give a good match of hybrid and baseline digital data along the sea-level, 
static operating line. For other flight conditions, the intercept of the compressor dis- 
charge specific-heat function and the time-invariant fan discharge (hub) specific heat 
were adjusted to match baseline rotor speeds at the military power setting. The follow- 
ing equations define the dynamics in the F100-PW-100 real-time simulation: 
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‘Adjusted for each flight condition to match rotor speeds at military power setting. 
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Adjusted for each flight condition to match rotor speeds at military power setting. 
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APPENDIX C 


DIGITAL PROGRAM 
FORTRAN Listing 


DIMENSION VALS(248),ZSC(2),PADR(3),PVAL(3),IDF0RM(4) 

C*****ADC VARIABLES FOR MAIN LOOP 

SCALED FRACTION X0 ,XI ,X2,X3,X4,X5,X6,X7,X8,X9,Xl 0,X1 1 
SCALED FRACTION XI2,Xl3,X!4,Xl5,Xl6,XlT,Xl8 
C*****ADC VARIABLES FOR OUTPUT LISTING 
SCALED FRACTION X19,X20,X2I 
C*****DAC VARIABLES 

SCALED FRACTION Y2,Y3,Y4,Y5, Y6,Y7,Y8,Y9,Y1 i,YI I ,Yl2 
C*****0THER VARIABLES 

SCALED FRACTION XVALSCI I ,8,6) , YVALS(8,6) ,ZVALS(1 0,8, 12) ,Yl (24) , 
1V9,V10,SSQRT,MAP2,MAP2L,V1 ,V2,V3,V4,V5,AR,PRSUB,DPR,PE, 
2PRE,PRSUP,FRD,V6,V7,V8,RT4,RT4I 
C OMMON/M APS /XVALS , YVALS , Z VALS , IX ( 6) , JY ( 6) , NX ( 6) , NY ( 6) , 
1KX(6),KY(6) 

LOGICAL SENSW 
REAL Ml 

C*****I NPUT REAL COMPONENT DATA AND MAP SCALE FACTORS 
TYPE I 
TYPE 101 

1 F0RMAT(/3X,33HPLACE DATA CARDS FOR MAPS NO. 2-9) 

101 FGRMAT(3X,34HAND DAC ICS IN READER. THEN R-S-R./) 

PAUSE 

C****4fiEAD COMPRESSOR SVA AND COMPRESSOR P/P MAPS 
DO 5 N=2,3 

READ (6*4) (IDFORM(I) ,1=1,4) ,NCV,NPT 
READ(6,2)XSC,YSC,ZSC(1) 

2 FORMAT (4F8.2) 

NTBL=NCV*(2*NPT+1) 

RE AD (6,3) (VALS (I ) , 1=1 ,NCV) 

3 FORMAT (8F8. 4) 

4 FORMAT MA2, 21 3) 

Jr NCV+1 

R EAD(6, IDFORM) (VALS ( I) , 1= J,NTBL) 

C*****SCALE MAPS 2-3 
DO 5 J=1,NCV 
YVALS (J , N)= VALS (J)/YSC 
DO 5 1=1, NPT 
JX=NCV+2*(J-1 )#NPT+I 
XVALS ( I , J , N ) = VALS ( JX ) /XSC 
JZ 1= JX+NPT 

ZVALS(I,J,N)=VALS(JZ1)/ZSC(1) 

IX(N)=1 
JY(N)= 1 
NX(N)=NPT 
NY(N)=NCV 
KX(N)r| 

KY(N)=0 

5 CONTINUE 

C*****flEAD FAN, HIGH TURBINE AND LOW TURBINE MAP DATA 
DO 7 N= 4 ,6 
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READ (6, 4) (IDFORM(I) ,1 = 1 ,4>,NCV,NPT 
READ(«,2)XSC,YSC,ZSC<1 ) ,ZSC(2) 
READ(6,3)<VALS(I),I=1,NCV> 

NTBL=NCV*(3*NPT+1 ) 

j - j 

R EAD <6,1 DFORM ) < VALS ( I ) , I = J , NTBL > 

C*****SCALE MAPS 4-9 
511 DO 6 J=1,NCV 

YVALS(J,N)=VALS(J)/YSC 
DO 6 1:1, NPT 
JX= NCV+3*(J- 1 )*NPT+I 
XVALS < I , J , N ) = VALS < JX ) /XSC 
JZUJX+NPT 

ZVALS (I , J, N>: VALS (JZl >/ZSC < 1 ) 

JZ2= JZ 1+NPT 
N2=N+3 

Z VALS (I , J , N2)= VALS < JZ2) /ZSC (2) 

6 CONTINUE 
IX<N>=1 
JYCN)=1 
NXC N>= NPT 
NY(N)= NCV 
KX(N):I 
KY<N):0 

7 CONTINUE 

C*****INITIALIZE dacs 

READ(6,9) <YI(I),I=1 ,24) 

9 FORMAT (<5<5X,S7>)) 

C*****INITIALIZE ANALOG CONSOLE 
CALL QSHYIN(IERR ,681,680) 

CALL QSRUN(IERR) 

CALL QSMSN(IERR) 

CALL QSC(I,IERR) 

CALL QSCd , IERR) 

DO If M= 1 ,24 
MM=5*M-1 

CALL QWFRL (MM,. FALSE., IERR) 

II CONTINUE 

C*****SPECIFY FLIGHT CONDITION FOR ENGINE 
C***** TRANSFER INITIAL DAC VALUES 
115 CALL QW8DAS (YI , I ,24 , IERR ) 

CALL QSTDA 
TYPE 11 
TYPE 111 

11 FORMAT </3X, 33HTYPE DESIRED VALUES FOR PI,TI,MI.) 
Ill FORMAT ( 3Xi 33HDEPRESS SSW(D) FOR NO INLET LOSS,/) 

ACCEPT 12,PI,T0 ,MI 

12 FOR MAT (F7 ,3,F8*3,F5.2) 

T2QI= (1,+.2*M0*M0) 

P2Q0: T2Q0**3.5 


ETAI= l.B 

IF(SENSW(4) > GO TO 13 

IFtMB.LE.l.) GO TO 13 

ETAI= 1.1 - . B 7 5* (M B- 1 . ) ** 1 . 35 

13 T2 = T2QB*T8 

P2r P2Qfl*PB * ETAI 
PBSI= PB*. 68948 
TBSI = TB*. 55555 
P2SI= P2*. 68948 
T2SI = T2*. 55555 
TYPE 14,P2,T2 

14 FORMAT C/3X»5HP2 = ,F7.3,6H T2 = ,Fs.3/) 

VI = P2/4B. 

V2 = T2/1BBB. 

V3= PB/2B. 

V5=MB/3. 

V6rTB/lBBB. 

V5=V5*SSQRT<V6> 

C*****PLACE ANALOG IN IC MODE 
TYPE 2 B 
TYPE 2B1 

2B FORMAT </3X,36HSLAVE CONSOLES 2 AND 3 TO CONSOLE 1.) 

2B 1 FORMAT <3X,18HMANUALLY GO TO IC.) 

TYPE 21 

21 FORMAT (3X,44HPR0CEED TO DYNAMIC PART OF PROGRAM BY R-S-R./) 
PAUSE 

C*****R EAD ADC VALUES AND GENERATE MAP OUTPUTS 

22 CALL QRBADS(XB,B,19'IERR) 

Y7=MAP2<4,XB,Xl> 

IF(X8.GT..BBBBBS) GO TO 221 
Y7=Y7*(.99999S+.B6633S*X8> 

221 Y5=MAP2L<7> 

1 F (X9 .GT..BBBBBS) X9=.BBSBBS 
V4=MAP2<2,X3,-X9) 

Y6=MAP2 <3,X2,X3) 

Y6=<.5B1BBS*Y6)/(.5BBBBS*V4+.5BBBBS) 

RT4-SSQRT <X17) 

V7=.86138S*X5/RT4 

Y8=MAP2(5,X4,V7) 

V9=MAP2L<8) 

RT41 = SSGRT(X18) 

V8=.84864S*X7/RT41 

Y1B=MAP2<6,X6,V8) 

V1B=MAP2L<9) 

Y9=-V9*RT4/.84895S 

Y11=-V1B*RT41/.5BBBBS 

Y2=Y6 

Y3=Y7 

Y4=X1B*X11 

Y5=Y5/.998BBS 
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AR=X12/.97838S-.6141IS*X12*X12+.26I76S 

DPR= < . 69683S*AR*AR- AR+.32I94S)/ . 1 I567S 

PRSUB=.52828S-DPR 

PE=PRSUB*Xl3 

FRD=.5I481S*X16*V5 

IF(.33333S*V3.LT»PE) QO TO 225 

PRE=.33333S*V3/X13 

Y12= (XI 5*SSQRT ( ( .2578S- .257BS*PRE)*Xl4>-FRD)/.34933S 
GO TO 23 

225 PRSUP= .52828S+DPR 

PE=PRSUP*X13/.66667S 

Y12= (X15*SSQRT< (.473I3S-PRSUP/.9I1 I6S+.B5I65S*PRSUP*PRSUP) 
1 *X14>-FRD-.931 54’S*X12*AR*( ,5IIIIS*V3-PE) )/.34933S 
C*****TRANSFER UPDATED DAC VALUES TO ANALOG 
23 CALL QVBDAS(Y2,2,11 ,IERR> 

CALL QSTDA 

C*****OUTPUT UNSCALED DATA AT TELETYPE IF DESIRED 
IF(.NOT.SENSV(l>> GO TO 22 
CALL QRBADS(X19,19,3,IERR> 

PI 3Q2 = XI 

P13Q2 s P13Q2*1 .5*3.1116 
XNLR2= XI 

XNLR2=XNLR2*1 .2*9872.8 
P3Q22=X2 

P3Q22= P3Q22*I. 5*8. 4621 
XNHR22 = X3 

XNHR22=XNHR22*1. 1*9942. 3 
P41Q4=X4 

P41Q4=P4 104*3. I*. 25784 
CNHPT = V7 

CNHPT=CNHPT*1 .25*221.28 
CNHTSI =CNHPT*1 .3416 
P5 04 1 =X6 

P5041=P5Q4l*2.l*.47763 
CNLPT = V8 

CNLPT=CNLPT*1 .5*215.14 
CNLTSI= CNLPT*1 .3416 
GVIP0S=X8 

GVIP0S=GVIP0S*25.- 5.4 
HVSP0S=X9 

HVSP0S=HVSP0S*4l.+2.5 

VF4=X1I 

VF4 = VF 4*4 . 

VF4SI= VF4*. 45359 
VF4=WF4*36II. 

AN-X12 

AN = AN*1 III. 

ANSI = AN*. 11164516 
ANr AN/144. 

P4= X2 1 


in 
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P4 = P4*5fl. 

P4SI = P4*. 68948 
XNL=X7 

XNL= XNL*1 51 If. 

XNH=X5 

XNH= XNH*15llf. 

T4=X17 

T4= T 4*411 I • 

T4SI = T4*. 55555 
PLA = X19 
PLAr PLA*15f • 

T7=XI4 

T7 = T7*5III. 

T7SI = T7*. 55555 
WF7 = X2f 
WF7 = WF7*2I. 

WF7SI= WF7*. 45359 
WF7=VF7*36fl. 

WA2=X16 
WA2= WA2*3fl. 

WA2SI= WA2*. 45359 
T41 = Xl8 
T41=T41*3III. 

T4iSI=T4l*. 55555 
P2 1 02= Y5 

P21Q2=P21Q2*1. 5*2. 9931 
WAR 22= Y6 

WAR22=WAR22*1 .1*54.988 
WR22S1= VAR22*. 45359 
WAR2 = Y7 

WAR2=WAR2*1 .2*221 .57*1,117 
WAR2SI= WAR2*. 45359 
WPHPT= Y8 

WPHPT=WPHPT*1 .31*. 18257 
WPHTSI= WPHPT*. 36548 
HPHPT = V9 

HPHPT=HPHPT*1 .25*. 31 143 
HPHTSIr HPHPT*3118.7 
WPLPT= Yll 

WPLPT=WPLPT*2 . 5* . 31 76 1 
WPLTSI= WPLPT*. 36548 
HPLPT : VII 

HPLPT= HPLPT*2 . 5* . 1 9522 
HPLTSI= HPLPT*31 18.7 
FN=Y12 
FN=FN*3lf If • 

FNSIrFN*4.4482E-3 
TYPE 24 

24 FORMAT (21X.33HF1 II SIMULATION STEADY-STATE DATA //> 
IF( , N0T.SENSWC2) ) GO TO 441 



TYPE 25,P0SI,P0 

25 FORMAT(5X»9HP0 = ,F7.3,9X,7HN/SQ CM,7X,2H< ,F7.3,9X,8HPSIA ) 

1 ) 

TYPE 26,T0SI,T0 

26 FORMAT(5x,9HT0 = ,F7.2,9X,7HK ,7X,2H( ,F7.2,9X,8HR ) 

1 ) 

TYPE 27, Mf 

27 FORMAT (5X.9HM0 = ,F7.4> 

TYPE 28,P2SI ,P2 

28 F0RMAT(5X,9HP2 = ,F7.3,9X,7HN/SQ CM,7X,2H< ,F7.3,9X,8HPSI A ) 
1 ) 

TYPE 29 T2SI T2 

29 FORMAT (5X»9HT2 = ,F7.2,9X,7HK ,7X,2H( ,F7.2,9X,8HR ) 

1 ) 

TYPE 3 I, PL A 

31 FORMAT (5X»9HPL A = ,F7.2,9X,3HDEG) 

TYPE 3I,VF4SI,WF4 

31^ FORMAT <5X»9HWF4 = ,F7.4,9X,7HKG/SEC ,7X,2H< ,F7. I,9X,8HLBM/HR ) 

TYPE 32, WF7SI, WF7 

32 FORMAT ( 5X»9HWF7 = ,F7.3,9X,7HKG/SEC ,7X,2M( ,F7. 0,9X,8HLBM/HR ) 
I) 

TYPE 33,GVIPOS 

33 FORMAT (5X»9HGVIP0S = ,F7.3,9X,3HDEG) 

TYPE 34 ,HVSPOS 

34 F0RMAT(5X,9HHVSP0S = ,F7.3,9X,3HDEG) 

TYPE 35, ANSI, AN 

35 FORMAT <5X,9HAN = ,F7 .5 ,9X,7HSQ M ,7X,2H< ,F7.4,9X,8HSQ FT ) 
1 ) 

TYPE 36,XNH 

36 F0RMAT(5X,9HXNH = ,F7.I,9X,3HRPM) 

TYPE 37 ,XNL 

37 FORMAT <5X,9HXNL = ,F7.0,9X,3HRPM> 

TYPE 38 WA2SI WA2 

38 F0RMAT(5X,9HWA2 = ,F7.2,9X,7HKG/SEC ,7X,2H( ,F7.2,9X,8HLBM/SEC> 

1 ) 

TYPE 40 ,P4SI , P4 

40 FORMAT (5X»9HP4 = ,F7.2,9X,7HN/SQ CM,7X,2H< ,F7.2,9X,8HPSIA ) 

1 ) 
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TYPE 41 , FNSI, FN 

41 F0RMAT(5X,9HFN = ,F7.2,9X,7HKN ,7X,2H< ,F7.0 ,9X,SHLBF ) 

1 ) 

TYPE 42,T4SI,T4 

42 FORMAT (5X»9HT4 = ,F7.1,9X,7HK ,7X,2HC ,F7.1,9X,8HR ) 

1) 

TYPE 39,T41SI,T41 

39 FORMAT <5X,9HT41 = ,F7.1,9X,7HK ,7X,2H< ,F7.1,9X,8HR ) 

1 ) 

TYPE 43,T7SI,T7 

43 FORMAT (5X»9HT7 = ,F7.1,9X,7HK ,7X,2H< ,F7.1,9X,8HR ) 

1 ) 

TYPE 45.P13Q2 

IF(.N0T.SENSW<3)> GO TO 22 


45 

FORMAT (5X»9HP13Q2 = 

,F7,4) 
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TYPE 46,XNLR2 





46 

FORMAT (5X»9HXNLR2 = 
TYPE 47,VTAR2SI , WAR2 

,F7.0,9X,3HRPM> 




47 

F0RMAT<5X,9HWAR2 = 

,F7.2,9X,7HK G/SEC 

,7X,2H< 

,F7 .2, 9X,8HLBM/SEC) 


1> 

TYPE 4B,P21Q2 
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FORMAT (5X.9HP21Q2 = 

TYPE 49,P3Q22 

,F7.4) 




49 

FORMAT <5X»9HP3Q22 = 

TYPE 50,XNHR22 

,F7. 3) 




50 

FORMAT <5X,9HXNHR22 = 
TYPE 51 , WR22SI ,WAR22 

,F7. 0 ,9X,3HRPM) 




51 

FORMAT <5X»9HWAR22 = 

,F7 ,2 ,9X, 7HKG/SEC 

,7X,2H( 

,F7,2,9X,8HLBM/SEC) 


1) 

TYPE 52,P41Q4 





52 

FORMAT (5X»9HP41Q4 = 

TYPE 53,CNHTSI,CNHPT 

,F7. 5) 




53 

FORMAT <5X,9HCNHPT = 

,F7.2,9X,7H 

,7X,2H( 

,F7.2,9X,8H 

) 


1) 

TYPE 54, WPHTSI ,WPHPT 





54 

FORMAT (5X,9HWPHPT = 

,F7.5,9X,7H 

,7X,2H( 

,F7.5,9X,8H 

) 


1) 

TYPE 55,HPHTSI ,HPHPT 





55 

FORMAT (5X,9HHPHPT = 

,F7 * 1 ,9X, 7H 

, 7X,2H< 

,F7.5,9X,8H 

) 


1) 

TYPE 56,P5Q41 





56 

FORMAT (5X,9HP5Q41 = 

,F7 ,5) 





TYPE 57,CNLTSI,CNLPT 
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57 FORMAT (5X»9HCNLPT 


,F7.2,9X,7H 


,7X,2H< ,F7.2,9X,8H 


) 


1 ) 

TYPE 58,WPLTSI,VPLPT 

58 FORMAT <5X*9HWPLPT s ,F7.5,9X,7H ,7X,2H< ,F7.5,9X,8H ) 

I) 

TYPE 59,HPLTSI,HPIPT 

59 FORMAT (5X.9HHPLPT = ,F7.J,9X,7H ,7X,2H< ,F7.5,9X,SH ) 

1) 

QO TO 22 
END 
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FORTRAN Symbols 


AN 

ANSI 

AR 

CNHPT 

CNHTSI 

CNLPT 

CNLTSI 

DPR 

ETAI 

FN 

FNSI 

FRD 

GVIPOS 

HPHPT 

HPHTSI 

HPLPT 

HPLTSI 

HVSPOS 

I 

ID FORM 
IERR 
IX 
J 

JX 

JY 

JZ1 

JZ2 

KX 


2 

exhaust nozzle throat area, ft 

2 

exhaust nozzle throat area, m 

exhaust nozzle expansion ratio (scaled) 

high -pressure -turbine corrected speed, rpm/°R^^ 

high -pressure -turbine corrected speed, rpm/K*^ 

low -pressure -turbine corrected speed, rpm/°R^^ 

low -pressure -turbine corrected speed, rpm/K^^ 

shift in critical pressure ratio due to expansion ratio 

inlet efficiency 

net thrust (uninstalled), lbf 

net thrust (uninstalled), kN 

ram drag (scaled) 

inlet guide vane position, deg 

high -pres sure -turbine enthalpy drop parameter, Btu/lbm-°R^^-rpm 
high -pressure -turbine enthalpy drop parameter, J/kg-K'^-rpm 
low -pres sure -turbine enthalpy drop parameter, Btu/lbm -°R -r pm 

low -pressure-turbine enthalpy drop parameter, J/kg-K^^-rpm 
stator vane position, deg 
integer 

array containing alphanumeric format data 
error flag for linkage routines 

array containing number of points per curve for each map pair 
integer 

map scaling index 

array containing number of curves for each map pair 
map scaling index 
map scaling index 

array containing x out -of -range counts for each map pair 



KY 

MAP 2 

MAP2L 

MM 

MO 

N 

NCV 

NPT 

NTBL 

NX 

NY 

N2 

PADR 

PE 

PLA 

PRE 

PRSUB 

PRSUP 

PVAL 

PI 

PL3I 

PJQI 

QRBADS 

QSC 

QSHYIN 

QSMSN 

QSPC 

QSRUN 

QSTDA 

QWBDAS 


array containing y out -of -range counts for each map pair 

bivariate function (first function) 

bivariate function (second function) 

function relay address 

Mach number 

map index 

number of curves for map being read in 

number of points per curve for map being read in 

number of data points for map being read in 

array containing number of points per curve for each map pair 

array containing number of curves for each map pair 

integer 

array containing potentiometer addresses 

exhaust plane pressure (scaled) 

power lever angle, deg 

nozzle pressure ratio 

critical nozzle pressure ratio 

design pressure ratio for supersonic nozzle flow 

array containing potentiometer settings 

pressure at station I, psia 

p 

pressure at station I, N/cm 

ratio of pressure at station J to pressure at station I 
linkage routine for reading ADC’s 
linkage routine for selecting analog console 
linkage routine for addressing analog consoles 

linkage routine for setting analog time scale to normal milliseconds (N-MS) 

linkage routine for selecting potentiometer coefficient mode 

linkage routine for selecting logic run mode 

linkage routine for transferring DAC data 

linkage routine for loading DAC 's 
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QWFRL linkage routine for positioning function relays 


RT4 square root of T4 (scaled) 

RT41 square root of T41 (scaled) 

SENSW array containing logical indication of sense switch positions 
SSQRT scaled -fraction square root routine 

TJ temperature at station I, °R 

TISI temperature of station I, K 

TJQI ratio of temperature at station J to temperature at station I 
VALS array containing unsealed map data 

VI fan inlet pressure (scaled) 

V2 fan inlet temperature (scaled) 

V3 ambient pressure (scaled) 

V4 shift in compressor map corrected airflow due to change in stator 

vane position (scaled) 

V5 Mach number (scaled) 

V6 ambient temperature (scaled) 

V7 high-pressure -turbine corrected speed (scaled) 

V8 low -pressure -turbine corrected speed (scaled) 

V9 high-pressure -turbine enthalpy drop parameter (scaled) 

V10 low -pres sure -turbine enthalpy drop parameter (scaled) 

WA2 fan airflow, lbm/sec 

WA2SI fan airflow, kg/sec 

WAR2 fan corrected airflow, lbm/sec 

WAR22 compressor corrected airflow, lbm/sec 

WAR2SI fan corrected airflow, kg/sec 

WF4 main -combustor fuel flow, lbm/hr 

WF4SI main -combustor fuel flow, kg/hr 

WF7 augmentor fuel flow, lbm/hr 

WF7SI augmentor fuel flow, kg /hr 

WPHPT high-pressure -turbine corrected flow, lbm -°R-in. ^/lbf-rpm -sec 
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WPHTSI 

WPLPT 

WPLTSI 

WR22SI 

XNH 

XNHR22 

XNL 

XNLR2 

XSC 

XVALS 

XI 

YSC 

YVALS 

YT 

ZSC 

ZVALS 


n 

high-pressure -turbine corrected flow, kg-K-cm /N-rpm-sec 

low -pressure-turbine corrected flow, lbm-°R-in. ^/lbf-rpm-sec 

o 

low -pres sure -turbine corrected flow, kg-K-cm /N-rpm-sec 

compressor corrected airflow, kg/sec 

high-speed-rotor speed, rpm 

compressor corrected speed, rpm 

low -speed -rotor speed, rpm 

fan corrected speed, rpm 

scale factor for map input x 

array containing scaled map input x data 

variable read on ADC channel I 

scale factor for map input y 

array containing scaled map input y data 

variable output of DAC channel I 

array containing scale factor for map output z 

array containing scaled map output z data 
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Compressor - Console 1 
















Hybrid interface - Console 1 







Control Interface - Console 1 

























Augmentor Volume - Console 2 
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Nozzle and Torque Balances - Console 2 














F100-PW-100 Real-Time Simulation Potentiometer Settings 
Used With 680 Analog Computer (Console 1)^ 


Poten- 

tiometer 

Setting 

Poten - 
tiometer 

Setting 


Poten- 

tiometer 

Setting 


Poten- 

tiometer 

Setting 

POO 

0.9330 

P30 

0. 5410 


P58 

0.5187 


P101 

0. 3993 

P01 

.6667 

P31 

.4000 


P59 

.8241 


P103 

. 1333 

P03 

.0001 

P33 

.4544 


P60 

.0768 


P105 

.0607 

P06 

.3006 

P35 

. 3876 


P61 

.0648 


P106 

. 1753 

P07 

.7274 

P36 

. 3109 


P62 

.0384 


P107 

. 1336 

P08 

.2000 

P37 

.7348 


P67 

.5000 


P108 

.0575 

P09 

. 8333 

P38 

. 1576 


P69 

.7642 


P110 

.2840 

Pll 

.2388 

P40 

. 89 85 


P70 

.3333 


Pill 

.2528 

P12 

.7430 

P42 

. 8333 


P73 

.8467 


P112 

.4996 

P13 

.0785 

P43 

.4109 


P75 

.3608 


P113 

. 8898 

P15 

.5000 

P45 

.6566 


P77 

.4699 


P115 

.9394 

P16 

. 1215 

P46 

.9891 


P84 

. 1333 


P117 

.0277 

P17 

.9010 

P47 

.7027 


P85 

.3296 


Q02 

.7603 

P18 

.2674 

P48 

. 5001 


P87 

. 1752 


Q04 

.0000 

P20 

. 1339 

P50 

. 1502 


P88 

.3001 


Q07 

.0000 

P21 

.2969 

P51 

. 3336 


P90 

. 1000 


Q09 

.9458 

P22 

.0606 

P52 

.0466 


P91 

.5780 


Q12 

.4867 

P23 

.4607 

P53 

.3674 


P93 

.4417 


Q14 

.4156 

P24 

.0712 

P55 

.4127 


P94 

.5011 


Q17 

.0767 

P25 

.3333 

P56 

.9132 


P95 

.2221 


Q19 

.6989 

P27 

.5442 

P57 

.4525 


P100 

.4845 


Q22 

.7823 

P28 

.4530 
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^Scale factors and parameters are in engineering units (U. S. customary units) . 



F100-PW-100 Real-Time Simulation Potentiometer Settings 
Used With 681 Analog Computer (Console 2)^ 


Poten- 

tiometer 

Setting 

Poten- 

tiometer 

Setting 

Poten- 

tiometer 

Setting 

Poten- 

tiometer 

Setting 

POO 

0. 8150 

P28 

0. 5000 

P79 

0. 8550 

P106 

0.2000 

P01 

.6667 

P29 

.2511 

P80 

.4414 

P107 

.2000 

P02 

. 1858 

P30 

.7897 

P81 

.5000 

P108 

.6667 

P03 

.2138 

P31 

.2789 

P82 

.7500 

P109 

.4300 

P04 

.0253 

P32 

.2000 

P83 

.3750 

P110 

.7054 

P05 

.4000 

P33 

.2000 

P84 

.2120 

Pill 

. 1874 

P06 

.7537 

P34 

. 5000 

P85 

.7470 

PI 12 

. 8950 

P07 

.5000 

P35 

. 6582 

P86 

.04 00 

P113 

.4140 

P08 

. 1000 

P36 

. 3333 

P87 

.0400 

P114 

.3361 

P09 

.2810 

P37 

. 3333 

P88 

. 1853 

P115 

.0023 

P10 

.6181 

P38 

. 5000 

P89 

.5000 

P116 

.6056 

Pll 

.7500 

P39 

.2687 

P90 

. 1415 

PI 17 

. 5144 

P12 

.2000 

P40 

. 1500 

P91 

.5517 

P118 

. 5113 

P13 

. 1332 

P41 

. 1500 

P92 

.1200 

P119 

.9999 

P14 

.7537 

P42 

.7897 

P93 

.4000 

Q02 

. 1885 

P15 

.8157 

P43 

. 6582 

P94 

.0640 

Q04 

.6212 

P16 

.7537 

P44 

. 5375 

P95 

.3358 

Q07 

.7957 

P17 

.3333 

P45 

. 6758 

P96 

.4280 

Q09 

.0000 

P18 

.3333 

P46 

. 6039 

P97 

.2608 

Q12 

.5384 

P19 

.2004 

P62 

.4000 

P98 

.0001 

Q14 

.3200 

P20 

.5562 

P63 

.0463 

P99 

.6667 

Q17 

.4197 

P21 

.8000 

P65 

.3333 

P100 

.0001 

Q19 

.6003 

P22 

.8000 

P67 

. 1486 

P101 

.8620 

Q22 

.0025 

P23 

.1753 

P72 

. 0634 

P102 

.3418 

Q24 

.3986 

P24 

.4664 

P73 

. 1336 

P103 

.6248 

Q27 

.6392 

P25 

.0001 

P77 

.0133 

P104 

.9999 

Q29 

. 1010 

P26 

.6667 

P78 

.2058 

P105 

.7080 



P27 

.3953 








Scale factors and parameters are in engineering units (U. S. customary units). 
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TABLE I. - DESIGN PARAMETERS 


Compressor discharge volume, V 3 , m 3 0.0468 

Main -combustor volume, V^, m 3 , 0 0468 

Interturbine volume, V 4 m3 0 6554 

Mixing volume, V g , m 3 . 0. 847 0 

Augmentor volume, V ? , m 3 0.7046 

Duct volume, V- 0 , m3 1 


Augmentor inductance, (Z/Ag) N-sec 2 /kg-cm 2 0 0007598 

^ AB 

Duct inductance, (J/Ag c )^, N-sec2/kg-cm2 0.007598 

High -speed -rotor inertia, I H , N-cm-sec 2 515.2 

Low-speed-rotor inertia, 1 ^, N-cm-sec 2 610.0 

Main-combustor pressure loss coefficient, Kg, N 2 -sec 2 /cm 4 -K-kg 2 0.00165 

Low-pressure -turbine discharge pressure loss coefficient, K pR5 1. 022 

Augmentor pressure loss coefficient K AB , N 2 -sec 2 /cm 4 -K-kg 2 4.3716xl0" 6 

Nozzle flow coefficient, K N , kg-K 1 // 2 /N-sec . 0. 1557 

Fraction of high-pressure -turbine cooling bleed that performs work, K ny lin _. 0.55 


_ BLWNT* 

Fraction of low-pressure -turbine cooling bleed that performs work, tttt rp ... 0. 06 

Fan inlet specific heat, c p 2 , JAg-K 1002 

Compressor specific heat,’c p 2 2 , JAg-K . 1003 

Compressor discharge specific heat, c p 3 , JAg-K 1043 

Main-combustor specific heat, c p 4 , J/kg-K ' , . , 1147 

Interturbine specific heat, c p 4 1 ’ JAg-K 1079 

Mixing -volume specific heat, ’c p 6 , JAg-K 1064 

Augmentor specific heat, c „, Vkg-K 


Duct inlet specific heat, c 


'» JAg-K 1004 


Duct discharge specific heat, c p lg , JAg-K 1030 

Main -combustor specific-heat ratio, y 4 a l 333 

Interturbine specific-heat ratio, y^ , 1.347 

Mixing -volume specific-heat ratio, Vg 1 381 

Augmentor specific -heat ratio, b l 332 

Compressor discharge temperature time constant, t 3 , sec 0. 05 

Duct temperature time constant, r 13 , sec 0.50 

Heating value of fuel, HVF, JAg 4.407xl0 7 

^Effectively decreased by a factor of 20 to match baseline digital data. 

Effectively decreased by a factor of 10 to increase simulation stability. 
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Figure 1. - Schematic representation of F100-PW-100 augmented turbofan engine. 
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lesign main 










Power lever angle, deg 

(e) Altitude, 22 555 m ; Mach number, 1. 5. (f) Altitude, 22 555 m : Mach number, 2 . 5. 

Figure 7. - Comparison of open-loop hybrid and baseline digital steady-state data for high -pres sure -turbine inlet temperature. 
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(b) Fan speed. 



.2 .3 ,4 .5 .6 .7 .8 .9 1 


Fraction of design main -com bus tor fuel flow. 

Ic) Main-combustor pressure. 

Figure IL - C f .iparison of closed-loop hybrid, baseline digital, and experimental steady-state data for 
nonaugmented mode. Sea-level, static, 297. 9 K day conditions. 
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Time, sec 

Figure 14 - Comparison of hybrid and baseline digital simulation transient data for idle-to-military power lever slam. Id le-area -reset on; sea-level, 
static conditions. 
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Figure 14 - Concluded 
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Figure 13. - Concluded. 















Fraction of design high -pressure-turbine Fraction of design high-pressure-turbine 



(b> Enthalpy drop parameter. 

Figure 19. - FIOO-PW-lOO high-pressure-turblne performance maps. 
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Normalized duct pressure drop, (P 





Figure 23. - FlOO-PW-lOO nozzle area reduction required to match baseline sea-level, static 
a ug mentor data. 
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NASA-Langley, 1975 
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